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ABSTRACT 

We propose a new hypothesis for the origin of protoplanetary discs with large inner 
holes (or gaps), so-called transition discs. Our gas disc model takes into account lay- 
ered accretion, in which poorly-ionized low-viscosity dead zones are sandwiched by 
high- viscosity surface layers, and photoevaporative winds induced by X-rays from the 
central stars. We find that a gap opens at a radius outside a dead zone, if the mass 
loss rate due to photoevaporative winds exceeds the mass accretion rate in the dead 
zone region. Since the dead zone survives even after the gap opens, mass accretion 
onto the central star continues for a long time. This model can reproduce large gap 
sizes and high mass accretion rates seen in observed transition discs. 
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1 INTRODUCTION 

Spatial distribution of dust in protoplanetary discs is mea- 
sured from disc continuum spectra in near-infrared to mil- 
limeter wavelengths whereas the gas density is measured di- 
rectly by line emission from such as CO and indirectly from 
Ha emission which is used for estimations of accretion rates 
of gas onto central stars. Infrared excess due to dust emis- 
sion disappears in 3 Myr f or roughly a half of discs an d for 
almost all discs in 10 Myr (|Haisch. Lada fc Ladall200lh. Gas 
in inn er discs disappears in similar time-scales ( Fedele et al.l 
but the life time of gas in outer discs is not clearly 

known. 

, Recent observations, such by the Spitzer infrared space- 
telescope, the Submillimeter Array, and the Subaru tele- 
scope, revealed discs in the dis p ersal phase, so-called 
transition discs (iMen'n et al.l I2OI0I : iThalmann et al.l I2OI0I : 
[Andrews et al.ll2011^ . They have optically thick outer discs 
whereas inner discs with sizes up to ~ 70 AU are opti- 
cally thin. This indicates depletion of dust in the inner 
discs. Some of transition discs show near-infrared excess in- 
dicating formation of gaps instead of holes and some others 
have weaker mid to far infrared emission in dicating optically 
thinner outer discs (|Muzerolle et al.ll2010l ). The fraction of 
transition discs in protoplanet ary discs increases with time 
jCurrie fc Sicilia-Aguii3l201ll '): 15-20 per cent at 1-2 Myr 
and more than 50 per cent at 5-8 Myr. A pu zzling thing 
is that a large fraction (75 per cent: lMen'n et al . (2010)) of 
transition discs exhibit gas accretion onto their central stars 
and their accretion rates are close to those for classical T 



Tauri stars, ~ 10"** M, 



(|Hartmann et al.lll99^ ). 
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There are three proposed explanations for the origin of 
transition discs but all of them do not seem to b e fully sat- 
isfactory dAndrews et al.ll201ll : IWilliams fc Cieza 2011). (1) 
Gap opening by photoevaporative winds (Gorti, DuUemond 
& HoUenbach 2009; Owen, Ercolano fc Clarke 2011). If the 
mass loss rate of photoevaporative winds exceeds the mass 
accretion rate of viscous evolution, a gap can open. How- 
ever, before a gap opens, the mass accretion rate largely 
decreases and becomes much less than 10~* Mq yr~^. This 
is the case even with X-ray photoev aporation with the mass 
loss rate as high as 10"* Mq yr"^ l|Owen et al.ll201ll ). The 
sizes of the inner holes due to photoevaporation cannot be 
larger than 20 AU when mass accretion onto the central star 
still remains. (2) Dust growth (DuUemond fc Dominik 2005; 
Tanaka, Himeno fc Ida 2005). This is expected particularly 
in inner discs as the models show that the growth time of 
dust has strong dependence on distance from the central 
star. However, dust growth alone cannot be responsible for 
transition discs, because the models for dust growth pre- 
dict smooth opacity changes with radius whereas observed 
transition discs show clear gaps of dust. (3) Gap opening 
by giant planets. IZhu et all (|201ll ) show that multiple giant 
planets can open a gap as wide as the largest holes observed. 
However, depletion of the surface density in the gap results 
in a large decrease of the mass accretion rate. 

The most probable source of the disc viscosity is the 
turbulent vi scosity caused by the ma gneto-rotational insta- 
bility (MRI: iBalbus fc HowLevI (119911 )). Since the ionization 
rate in the thick inner disc is not high enough, a poorly ion- 
ized layer with a low viscosity called a dead zone forms near 
the disc midplane sandwiched by MRI active layers (Gam- 
mie 1996; Zhu, Hartmann fc Gammie 2010). This type of 
accretion is named layered accretion. In this paper, we de- 
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velop a gas disc model which takes into account layered ac- 
cretion and photoevaporative winds. Layered accretion has 
not been taken into account in previous works of photoe- 
vaporation. As we will show, our disc model can reproduce 
some of important properties of observed transition discs, 
such as large gap sizes and high mass accretion rates. 

In Sec. 2, we introduce our gas disc model. In Sec. 3, 
results are shown for both discs with and without layered 
accretion over wide ranges of parameters. In Sec. 4, implica- 
tions for transitions discs are discussed. In Sec. 5, the sum- 
mary of this study is given. 



2 METHODS 

2.1 Basic equations 

We assume that mass accretion from a molecular cloud core 
has already completed and that the mass of the central star 
is fixed during disc evolution. Time evolution of the surface 
density E is derived from the mass and angular momentum 
conservation equations; 



9E _ J_dM_ _ 
dt 2'Kr dr 

dm 1 d 



dt 



2-Kr dr 



Mj + Jvis) 



JwJ, 



(1) 

(2) 



where M is the mass accretion rate (assumed to be positive 
for inward accretion) , Sw is the mass loss flux of photoevapo- 
rative winds, j — r^Q is the specific angular momentum with 
Q being the Keplerian frequency, and Jvis = 2-Kr'^Tivd£l/ dr 
is the angular momentum flux transferred by the disc vis- 
cosity, v. From Eqs. ([l]) and ([2]), the mass accretion rate is 
obtained as 
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'OT^r —(l^vr ) 
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where Vr is the radial velocity. Substituting this form 
to Eg. ([Tjl. we obtain the w ell known diffusion equation 
l|Lynden^ell fc Pringle|[l973 ): 
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(4) 



2.2 Disc viscosity 

Using the viscosity parameter a (jShakura fc SunvaevlllQT^ ) , 
the viscosity is given as 

aRT 



ach 



(5) 



where c is the isothermal sound velocity at the disc mid- 
plane and h = c/JIk is the vertical scale height of the disc. 
The isothermal sound velocity is given as c = {RT/ fj,)^'''^ , 
where R is the gas constant, T is the disc temperature at 
the midplane, and we take the mean molecular weight as /i 
= 2.33 g moP^. For models without dead zones, a is radi- 
ally constant. For models with dead zones, the a parameters 
for active surface layers and dead zones are defined to be 
and Od. Then, the effective disc viscosity parameter is given 
as 



(6) 



Table 1. Rosseland mean opacity from [Stepinskil l|l99^ 



Applicability 


KO (cm^ g 1) 


/3 


T < 150 K 


2 X 10"** 


2 


150 K < r < 180 K 


1.15 X IQi* 


-8 


180 K < T < 1380 K 


2.13 X 10-2 


3/4 


1380 K < T 


4.38 X 10-** 


-14 



where Ea is the sum of the surface densities of both sides 
of active layers. In simulations, we adopt Oa = 0.01. The 
surface density Ea is modeled after Suzuki, Muto fc Inutsuka 
(2010): 



Ea — min 



EcR -I- Ex 



(lAu) 



(7) 



where Ecr is the surface density of layers ionized by cosmic 
rays and Ex is the surface density of layers ionized by X- 
rays from the central star at 1 AU. Equations (|6]) and ((Tjl 
indicate that a <C Qa if E ^ Ea. On the other hand, if 
E < EcR-l-Ex(r/lAU)"^, the region is entirely MRI active 
(E = Ea) so that a = eta. 

The surface densities Ecr and Ex particularly depend 
on dust amount in the disc, as reco mbination of ions and 
electrons occurs on the dust surfac e dSano et al ] |200Gl : [Bail 
[ioil). For example, ISuzuki et af] l|2010l) uses Eqr = 12 
g cm-'^ and Ex = 25 g cm-'^, whereas iGammij (Il996l ) 
adopts nominal values as Ecr = 200 g cm ~ ^ and Ex ~ 
0. The value of Ecr adopted in iGammij (|l996t ) comes 
from the attenuation l e ngth o f comic rays (~ 96 g cm-^; 
lUmcbavas hi fc Nakanj (j 198 il l). If the dust abundance is 
negligible, MRI is sustained in the layers as thick as the 
attenuation length (|Sano et al.ll200Gl). On the other hand, 
in the MHD simulations of lSuzuki et al.l ||2010|) . they a dopt 
the ionization degree calculated in ISano et all l|2000l ) and 
llnutsuka fc Sand (|2005l ) assuming that the dust-to-gas ratio 
be 0.01. In the present study, we vary Ecr as a parameter. 
On the other hand, we fix Ex to be 25 g cm^^, as we find 
that Ex is less important than Ecr for overall disk evolu- 
tion, in particular, for ad ^ 10~^ . MHD simulations show 
that dead zones have no n-zero viscosities l|Fleming fc Ston3 
[2003; ISuzuki et al.l l2O10l): the value of Qd is on the order of 
~ 10"^-10"^ 



2.3 Disc temperature 

As heat sources of discs, we take into account viscous heat- 
ing, irradiation from the central star, and background ir- 
radiation. We simply superpose all contributions which are 
independently modeled: 



: I irr ' amb ! 



(8) 



where Tvis, Th-y, and Tamb are the temperatures contributed 
from viscous heating, irradiation from the central star, and 
background irradiation. If only viscous heating is considered, 
the radiative cooling from both sides of the disc balances 
with the viscous heating as 



2<Ts: 



(9) 
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where Tvis^efF is the effective disc surface temperature. The 
refation between the midpfane temperature Ty ia and the ef- 
fective temperature is given bv lHubenvl (|l99(]l l as 



-'vis '^cff J vis.cff 7 
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3r 



1 + 



V3 
4 



1 

37 



(10) 



(11) 



where r = kE/2 is the optical depth at the midplane. 
Eq. pip is apphcable to all cases. For pure MRI active re- 
gions without dead zones, we set E = Sa and a = Qa. For 
the dead zone region with ad = 0, we set S q = EaQa s o that 
Eq. dTDI becomes the same with Eq. (7) of lGammi"^ l{199(i ) 
for rEa/E 3> 1. A detailed derivation of the first term in 
th e square bracket of the right hand side of Eq. (|ll|l is given 
in IWiinsch et al. I ||200Q1 . The opacity k is repre sented by 
K = k qT^ , where kq is a constant. We use k from IStepinskil 
shown in Table 1. Resuhs of the present study do not 
sensitively depend on opacity law. 

The disc eff ective temperature irrad iated by the central 
star is given by (|Ruden fc PoIlack|[l99ll ) 



3^ \V) ^ 2 \V) \r) [ 



din ft 
dinr 



- 1 



,(12) 



where T, and r* are the photospheric temperature and ra- 
dius of the central star. We adopt T, — 4000 K and r, — 3 
vq. The scale height h used in Eq. (|12p is the photospheric 
disc scale height and we simply assume it to be the same with 
the pressure scale height. After iHueso fc Guillod (|2005l ). we 
assume dlnh/dlnr = 9/7. This is validated as T oc r~^/^ at 
large radii where irradiation from the central star is usually 
the dominant heat source. The midplane temperature, Tin, 
irradiated by the central star is ca lculated using the vertica l 
temperature structure derived bv lMalbet fc BertoutI (|l99ll ) 
(see also Malbet, Lachaume & Monin (2001)): 
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T 
MO 



1 / -r 
+ - — exp 

4^0 V MO 



(13) 



Here mo is the averaged cosine of the incident angle of the 
central star given as 



Mo = 



JO ' 



(14) 



with the zeroth and first order moments, j" and of the 
incident intensity at the upper most layer: 

l/2l 

1-(1- 



J' 



47r 



47r 



(15) 



(16) 



Eq. (|13|l represents that Ti^r of an optically thick disc is 
lower than Ti^r.eff by a factor of 2^^^, as the upper most lay- 
ers emit half of the received energy toward the midplane. 
In the optically thin limit (r — >■ 0), the third term in the 
bracket dominates and we obtain T^^ = Tt{r,/rf/8, which 
corresponds to the equilibrium temperature irradiated by 



the half hemisphere of the central star without any obsta- 
cle. The temperature at r — > should be higher by a factor 
of 2^/'', but this difference does not affect our results as the 
stage with a low surface density is very short with photoe- 
vaporative winds. 

The ambient temperature is simply assumed as Tamb = 
20 K. 



2.4 Mass loss due to photoevaporation 

Among various types of photoevaporations, the one induced 
by X-rays from the central star is most strong and most 
import ant for gas disc evolution in th e planet formation 
regio n (|Ercolano. Clarke fc Drak^ l2009l : lOwen et all |2010| . 
l201ll). unless discs are close to very massive e xternal stars 
(| Adams et al] |2004| : iMitchell fc Stewart! 1201 d ). The total 
mass loss rate from a disc due to X-ray photoevaporation 
is estimated bv lOwen et all l|201ll ) as 



. = 6.4x10-«a(— ^ 

\ 10'='" ergs^ 



1.14 



Mq yr \ 



(17) 



where ^ is a constant of the order of unity (we fix it to 
be unity) and Lx is the X-ray luminosity from the cen- 
tral star. For sub-solar to solar m ass protostars, Lx is 10^®- 
10^^ erg s~MGudel et al.|[2007l) . Numerical simulations of 
lOwen et al.l (|2010l ) show that Ew is roughly proportional to 
^-3/2 ^^jjgjj. pjg Thus, we give the mass loss rate per 
unit area as 



-3/2 



47r(; 



.1/2 



1/2 



(for rx,in < J- < rx.out), (18) 



where rx.in and rx,out are the inner and outer edges of the 
range where mass loss occurs. We t ake r x.in ~ 1 AU and 
rx.out = 70 AU from I Owen et al] (|2010l ). Outside '^X.out , 
other mechanisms such as EUV or FUV photoevaporation 
l|Gorti et al.ll2009l ) may work, although these effects are not 
taken into account in the present work. Results of disc evo- 
lution do not sensitively depend on r-dependence of Ew, as 
similar results are obtained even if we vary the power-law 
index between -1 and -2. 

In some test simula tions, we also adop t the mass loss 
due to MRI disc winds (jSuzuki et al.ll2010l ). The mass loss 
rate of MRI disc winds at a given radius is proportional to 
the surface density. We find that gas depletes very rapidly 
even near the inner edge, and that the mass accretion rate 
onto the central star b ecomes less th a n lO' ^Mp yr"'^ within 
1 Myr, as shown in ISuzuki et al.l (|2010l ). This accretion 
rate looks too low as compared with t hose for typical discs 
around classical T Tauri stars ( Hartina nn et al.l|l9 98l . Prob- 
ably, the mass loss rate of MRI disc winds is much smaller, 
but this mechanism may still play an important role in disc 
evolution. 



2.5 Numerical procedures 

We numerically solve Eq. Q using the method described 
in Bath and Pringle (1981), in which the radial grid size 
is proportional to y'r. We use 1000 grids between 0.03 AU 
and 3000 AU. These numbers give the inner most grid size 
of 0.02 AU. At the inner boundary E is fixed to be zero, 
and at the outer boundary we adopt the outward mass fiux 
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Table 2. Input parameters for simulations of discs without dead 
zones. 



Run ID 


(1030 erg s-i) 


a 


Nl 


1.0 


10-2 


N2 




10-3 


N3 




3 X 10-" 


N4 




10-4 


N5 


3.0 


10-2 


N6 




10-3 


N7 




3 X 10-4 


N8 




10-4 



Table 3. Input parameters for simulations of discs with dead 
zones. We adopt «a = 0.01. 



Run ID Lx ScR ad 

(1030 erg s-^) (g cm-2) 



Dl 


1.0 


50 





D2 




100 




D3 


3.0 


50 




D4 




100 




D5 


1.0 


25 


10-5 


D6 




50 




D7 




100 




D8 


3.0 


25 




D9 




50 




DIO 




100 


); 


Dll 




200 




D12 


0.3 


25 


)j 


D13 




50 




D14 




100 




D15 


1.0 


50 


10-4 


D16 


3.0 






D17 


1.0 




10-3 


D18 


3.0 







given by —Stt'Ev. We find tliat the evolution of the total 
mass is nearly the same even with larger outer boundary 
radii whereas artificial mass loss from the outer boundary is 
not negligible with smaller outer boundary radii. The tem- 
perature is calculated by the Newton-Raphson method. 

For all simulations, the mass of the central star is the 
solar mass, Mq, the initial mass of the disc is 0.1 Mq, and 
the initial surface density is proportional to with the 
outer edge of 20 AU. The adopted time step size is 1 yr. A 
single simulation for ~ 10 Myr takes roughly one cpu day. 
We conduct 8 runs for discs with radially constant q's and 
18 runs for discs with dead zones. The input parameters are 
shown in Tables 2 and 3 for simulations of discs without and 
with dead zones. 



3 RESULTS 

Results for all simulations are summarized in Table 4. 
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Figure 1. Evolution of surface density, temperature, and mass 
accretion rate for run N2. The viscous a is 10-3 and is radially 
constant. Black dotted lines are initial conditions. Red, orange, 
blue, black, and purple solid lines are values at 0.1, 1.0, 5.0, 6.01, 
and 7.0 Myr from the beginning, respectively. Mass accretion onto 
the central star stops when the inner disc disappears at 6.12 Myr. 
In the panel of M, dashed lines represent for outward motion of 
gas while solid lines for inward motion. 



3.1 Discs without dead zones 

Figure 1 shows time evolution of the surface density E, the 
temperature T, and the mass accretion rate M for a disc 
with a radially constant a (run N2). As the disc spreads 
with time, E and T decrease. The kinks seen in the T profile 
correspond to the opacity transition temperatures (150 K, 
180 K, and 1380 K) and the corresponding kinks are also 
seen in E. The direction of the radial gas motion is inward 
in the inner disc and outward in the outer disc, and the 
zero radial velocity radius expands with disc expansion. The 
accretion rate M is independent of r in the inner disc and is 
equivalent to the mass accretion rate onto the central star 
M*, as long as M is sufficiently larger than the mass loss 
rate due to photoevaporation Mw When M^, becomes less 
than Mw, E rapidly decreases and eventually a gap opens at 
slightly beyond 1 AU, as well as simulations in lOwen et al.l 
(2011). Once a gap opens, the inner disc (inside the gap) 
disperses very quickly and mass accretion onto the central 
star stops. The hole size is only 1.3 AU at this time. The 
temperature increases when the disc becomes optically thin, 
but this phase (with a non-zero gas density) is very short. 

The gap opening slightly outside rx,in is explained as 
follows. Let us assume that the inward mass flux due to vis- 
cous accretion at rx.out is Mout- As gas moves inward by 
dr, the mass flux is reduced by 2irrT,vdr. Therefore, the in- 
ward mass flux decreases with decreasing r and eventually 
becomes zero at rx.in if Mout = M^. Thus, if M is indepen- 
dent of r without photoevaporation, a gap inevitably opens 
near rx.in- Since it takes a time for gas to radially move from 
»'x,out to rx,in (the viscous time scale at rx.out), Mout be- 
comes much less than Mw at the time of actual gap opening 
(see the black curve in Fig. 1). 

In some runs, large gaps can open (runs N3 and N6; 
see Table 4). This happens because the disc sizes are not 
sufficiently larger than rx,out and the mass accretion rates 
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Table 4. Summary of disc properties. 



Run ID 


A^*,lMyr 






^gap 


''"gap 


Mgup 


-^-^gap,out 


-^*,gap 




(IO-SMq yr-1) 


(Myr) 


(Myr) 


(Myr) 


(AU) 




(IO-^Mq) 


(10-9Mq yr-1) 


Nl 


7.73 


2.40 


- 


2.40 


1.47 


10.23 


10.23 


5.27 xlO"* 


N2 


9.19 


6.12 


- 


6.01 


1.34 


10.72 


10.72 


2.58 xlO-2 


N3 


6.73 


10.21 


- 


9.08 


60.37 


6.85 


1.25 


0.17 


N4 


3.22 


16.02 


- 


14.90 


1.34 


2.65 


2.61 


3.21 xlO-2 


N5 


1.27 


1.11 


- 


1.09 


1.22 


12.63 


12.62 


5.86 xlO-2 


N6 


6.80 


2.84 


- 


2.63 


60.37 


5.46 


1.37 


0.20 


N7 


5.31 


4.58 


- 


4.12 


1.34 


5.62 


5.59 


7.99 xlO-2 


N8 


2.56 


6.46 


- 


5.12 


1.34 


14.29 


14.21 


7.12 xlO-2 


Dl 


2.89 


12.17 


12.11 


4.36 


47.55 


47.37 


2.87 


2.89 


D2 


O.Ul 


o.u± 


7 Q9 


71 


'^l 1 s 


'^A 7Q 


A 77 


^ ^1 


D3 


2 89 


6.58 


6.54 


1.48 


54.58 


57.09 


2.12 


2.89 


D4 


5.61 


5.08 


5.03 


1.66 


42.42 


44.84 


3.94 


5.51 


D5 


3.48 


9.39 


9.39 


5.52 


57.02 


34.43 


1.85 


4.40 


UD 


D.oy 




O.lo 


c; Q7 


40. OZ 


1 n OS 
iU.yfs 


Q 1 7 
O.LI 


K fin 


D7 


18.31 


3.69 


3.46 


3.68 


2.05 


7.11 


7.11 


1.74 xlO-3 


D8 


3.26 


5.40 


5.37 


1.30 


61.23 


63.56 


1.03 


3.22 


D9 


6.53 


3.97 


3.94 


1.62 


57.02 


48.91 


2.10 


7.28 


DIO 


17.71 


2.51 


2.47 


1.83 


45.32 


20.08 


4.10 


15.88 


Dll 


30.85 


1.50 


1.37 


1.49 


2.54 


8.19 


8.19 


1.89 xlO-2 


D12 


3.55 


14.32 


14.15 


14.29 


28.20 


1.66 


1.65 


4.51 xlO-2 


D13 


7.00 


9.01 


8.28 


9.00 


1.22 


4.32 


4.32 


2.24 xlO-3 


D14 


18.49 


6.50 


4.07 


6.49 


1.10 


5.75 


5.75 


1.69 xlO-3 


D15 


17.37 


4.10 


3.95 


4.09 


2.05 


6.40 


6.40 


5.77 xlO-3 


D16 


15.84 


2.62 


2.58 


2.02 


52.18 


14.87 


2.98 


10.00 


D17 


19.99 


2.83 


2.20 


2.83 


1.47 


9.51 


9.51 


8.74 xlO"* 


D18 


15.31 


1.56 


1.48 


1.55 


2.72 


8.17 


8.17 


1.23 xlQ-^ 



^*,lMyr is the mass accretion rate onto the central star at i = 1 Myr, is the time when M« becomes 10~^^ Mq yr~^, is the hfe 
time of the dead zone (when a disc becomes entirely MRI active for all radii), igap is the time when a gap opens, r-gap is the gap radius 
where S becomes zero in the earliest time, A/gap and A/gap, out are the total disc mass and the disc mass outside rgap at t = tgap, and 

A/,, gap is Af* at t = tgap. 



are not radially constant near the outer edges of the discs. 
Even in these cases, the mass accretion rates onto the central 
stars at the time of gap opening are much smaller than those 
seen in transition discs. Overall, a large accretion rate and 
a large hole size are not simultaneously reproduced as long 
as a is radially constant. 

3.2 Discs with dead zones 

3.2.1 Cases with no dead zone viscosity (a^ = 0) 

Pictures of evolution of discs with dead zones are very dif- 
ferent from those without dead zones. Figure 2 shows evo- 
lution of a layered disc with Od ~ (run Dl). The mass 
accretion rate discontinuously drops at T = 150 K, be- 
cause the opacity law changes. As a result, two dead zone 
regions appear. Splittin g int o multiple dead zones was theo- 
retically predicted in ,Gammiei (1998) . Most of the disc mass 
remains in the dead zones and Aif* remains nearly constant 
as long as the dead zones exist. Since M in th e dead zone re- 
gion increases w ith r, E increases with time l|Gammiell 19961 : 
IZhu et al-llioiol ). 

The zero radial velocity radius initially locates at the 
outer edge of the outer dead zone. The direction of radial 
motion of gas is outward outside the outer dead zone, while it 
is inward in the dead zone region. The absolute magnitudes 
of the mass fluxes in both directions are similar to each 




0-1 1.0 10,0 100.0 1000.0 

Radius (AU) 



Figure 2. Same with Fig. 1 but for run Dl. A layered accretion 
is taken into account with = 0. Black dotted lines are initial 
conditions. Red, orange, blue, black solid lines are values at 0.1, 
2.0, 8.0, and 12.0 Myr from the beginning, respectively. The inner 
disc completely disappears at 12.17 Myr. 



other. Since the size of the dead zone is ~ 10 AU, most 
of gas removal by photoevaporation is taken place outside 
the dead zone region. A discussion of gap opening due to 
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Figure 3. Same with Fig. 1 but for run D6. A layered accretion is 
taken into account with = 10~^. Black dotted lines are initial 
conditions. Rod, orange, blue, black solid lines are values at 0.1, 
2.0, 6.0, and 6.22 Myr from the beginning, respectively. The inner 
disc completely disappears at 6.24 Myr. 



photoevaporation is similar to the case of Fig. 1, but now 
the mass flux is outward. Thus, a gap opens shghtly inside 
rx,out once M near the outer edge of the dead zone becomes 
less than Afw 

For run Dl, a gap opens at 48 AU roughly when the 
outer dead zone disappears, because M of the outer dead is 
larger than M„ while M in the inner dead zone is compa- 
rable to Mw. Since the inner dead zone survives for a long 
time, M« remains high even after a gap opens. The outward 
mass flux from the dead zone remains to be ~ Mw after a 
gap opens, as long as the dead zone exists. Once the dead 
zone disappears, the inner disc quickly dissipates and mass 
accretion onto the central star is halted. 



3.2.2 Cases with finite dead zone viscosity (a^ > 0) 

It is possible to retain a steady state mass accretion with a 
finite residual viscosity in the dead zone. The condition to 
retain a radially constant M is given from Eqs Q, (Oi a^nd 
(m as 



1/2 



d_ 

dr 



[[EaOa + (S - Ea)ad]rr^] = const. 



(19) 



If the radial profile of T is given such as irradiated discs, 
the above condition is fulfilled when Eod > SaCta so that 
the radial profile of S is adjusted. On the other hand, if 
viscous heating is the dominant heat source, the midplane 
temperature is given from Eqs. ((5]), (O, (|9]), (fTO|l as 



T = 



27 Rkq 
128 asBjJ, 



1 1/(3-/3) 



Sa)ad -f SaOa 



(20) 



^ Without photoevaporation, the direction of radial motion of 
gas outside the dead zone is outward onl y in the early ex pansion 
phase, but eventually turns to be inward l lZhu et al .l201(t) . except 
near the outer edge of the disc. The outward flux is retained by 
strong photoevaporation in our simulations. 



Figure 4. Radial locations of gaps. Left: rgap vs. for Scr = 
50 g cm~i (runs Dl, D3, D6, D9, and D15-D18). Asterisks and 
triangles are for cases of Lx = 1-0 and 3.0 in units of 10'^'' erg 
s~^. The values of rgap for Od = are shown at ad = 10~®. Left: 
rgap vs. ScR for = 10~^ (runs D5-D14). Diamonds, asterisks, 
and triangles are for cases of Lx = 0.3, 1.0, and 3.0 in units of 



10'^ 



respectively. 



Thus, if S^Qd }t SaQa, the radial profiles of T and E are 
mutually adjusted so that M can be independent of r. 

Figure 3 shows an example of evolution of a disc with 
Qd ~ 10~^ (run D6). Since this case satisfies the condition 
M becomes almost independent of r (see the 
orange line at 2.0 Myr in Fig. 3). Split into multiple dead 
zones is also suppressed. As well as the cases with Ud ~ 0, 
a gap opens at a radius outside the dead zone. 

Figure 4 shows the radial location of a gap rgap for 
various values of ad, Ecr, and Lx- We define rgap as a radius 
where E becomes zero in the earliest time while E > in 
outer radii. If A4 is independent of r in the dead zone region, 
the condition for gap opening is given by Afw A'/*- For 
small values of Qd and Ecr, Af, can be smaller than Afw 
even when a dead zone exists. Thus, a gap opens beyond 
a dead zone in such a case. On the other hand, for a case 
with large ad and Ecr, a gap opens at a small radius only 
after a dead zone disappears and M^, becomes sufficiently 
small, as well as the cases with radially constant q's shown 
in Section 3.1. Not surprisingly, gap opening beyond dead 
zones is possible even with large values of ad and Ecr if Afw 
(or Lx) is large. 



4 DISCUSSION 

Figure 5 shows comparison between observed transition 
discs and discs from our simulations. In the upper panel 
of Fig. 5, gap sizes and mass accretion rates are plotted. For 
simulations, we use rgap and the mass accretion rate onto the 
central star Af,,gap at the time of gap opening (see Table 4). 
Subsequent evolution curves are also shown. As can be seen, 
large gap (or hole) sizes and large accretion rates seen in ob- 
served transition discs are well reproduced in models with 
dead zones. On the other hand, models without dead zones 
are able to reproduce only discs with small mass accretion 
rates and small gap sizes. In runs N3 and N6, a gap opens 
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Figure 5. Comparison between modeled discs and observed tran- 
sition discs. Top: Gap (or liole) size vs. mass accretion rate onto 
the central star. Red squares represent observed transition discs, 
black and blue asterisks are all simulation results for discs without 
dead zones and with dead zones. Subsequent evolution curves are 
also shown. Small crosses represent subsequent evolution of Tgap 
and M« with an interval of 10^ yr only for runs N2, D l, and D6. 
Obser v ed data are f r om lEspaillat et al.l l20oi. l201(t) iKim et al.l 
1I2OO9I '). iMerm et akl ||2010|' |. and [Andrews et al.l | |2011^ . Bottom: 
disc mass vs. mass accretion rate onto the central star. Sym- 
bols are the same with those in the top panel. For simulations, 
we take the disc mass Mgap.out outs ide rgap, not the total disc 
mass, Mgap. Observational data from lNaiita. Strom Sz Muzeroll3 
are also added. 



at ~ 60 AU, but the mass accretion rate is small (in these 
runs, another gap opens at ~ lAU, and the disk outside the 
outer gap quickly dissipates). 

In the lower panel of Fig. 5, we compare between disc 
masses from observations and our simulations. For our sim- 
ulations, we take the disc mass outside the gap (Mgap, out 
in Table 4), not the total disc mass (Mgap), for consistency 
with observations. It should be noted that submillimeter - 
millimeter observations measure emission from dust primar- 
ily around ~ 100 AU and the surface density of gas is esti- 
mated using the interstellar gas-to-dust ratio. We find that 
our model cannot reproduce the masses of the most massive 
transition discs with high accretion rates, although the total 
disc masses including the inner discs in our simulations are 
comparable to the observed values (see Table 4). This means 



^ before a gap opens but it decreases 
after a gap opens (see Figs. 2 and 



Figure 6. Expected fraction of gapped discs for the cases without 
(left) and with (right) dead zones. Only discs with mass accretion 
onto their central stars are considered. The time-scales and tgap 
are the times when accretion onto the central star stops and when 
a gap opens. The luminosity Lx is 0.3 (diamonds), 1.0 (asterisks), 
and 3.0 (triangles) in units of 10"^" erg s~^. For runs with dead 
zone, Q(j = 10~^ (runs D5-D14). 



that most of the mass is in the dead zone for a disc with lay- 
ered accretion. In our model of layered accretion discs, E at 
100 AU is ~ 1 g cm 
to less than ~ 0.1 g cm 
3). On the other hand, observed transition discs sho w that 
E ~ 1-10 g cm"^ at 100 AU ([Andrews et al.ll201ll '). This 
value is similar to or even larger than those seen in classica l 
T Tauri discs (|Kitamura et al.||2002l : [Andrews et al.ll2010l '). 
The apparent contradiction between modeled and observed 
disc masses might be resolved if dust is accumulated near 
the inner edge of the outer d isc after a gap opens, resultin g 
in a large dust-to-gas ratio (| Alexander fc Armitagd [20071 '). 
The dust-to-gas ratio may also increase if photoevaporative 
winds remove only gas but not dust. 

The expected fraction of gapped discs is much larger 
for discs with dead zones than that for discs without dead 
zones (Fig. 6) , if the disc life times are similar. The fraction 
of transition discs in protoplanetary discs is more than 50 
per cent at 5-8 Myr l|Currie fc Sicilia-Aguil ar 2011), and a 
large fraction of them are accreting (jMerm et al..,2010 ). The 
fraction of transition discs with large holes (> 15 AU) is at 
least 20 per cent am ong millimeter-bright disc population 
(An drews et al.|[201ll ) although discs wi th large holes seem 
to be rarer than those with small holes ([Meri'n et al.ll20l"ol ). 
For discs with dead zones, it is possible to reproduce a high 
fraction of discs with large holes, if Eqr and are small and 
Lx is large. A detailed comparison may be able to constrain 
EcR and with a given distribution of Lx- 

Dead zones in our model are optically thick even after 
gaps open in discs. This is inconsistent with observations, 
although some of transition discs have inner optically thick 
region s indicated from near -infrared excess (pre-transition 
discs; lEspaiUat et~all (I2OIOI ')'). This indicates that dust re- 
moval due to their growth or migration is necessary whereas 
our current model adopts a fixed opacity. Since dust is usu- 
ally supplied from the outer discs (Brauer, DuUemond & 
Henning 2008), gap opening may make the inner discs or 
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dead zones optically thin. Alte rnatively, giant p l anets may 
clean up dust in the inner discs l|Merfn et ahlbOldlZhu et al.l 
Giant planets may also be responsible for transition 
discs with small holes and large accretion rates, as our model 
does not reproduce such discs (Fig. 5). 



5 SUMMARY 

In this study, we developed a gas disc model which takes 
into account layered accretion and photoevaporative winds 
induced by X-rays from the central stars. We found that a 
gap opens at a radius outside a poorly-ionized dead zone, 
if the mass loss rate due to photoevaporation exceeds the 
mass accretion rate in the dead zone region. Since the dead 
zone survives even after the gap opens, high mass accretion 
onto the central star remains for a long time. 

We found good agreements between modeled and ob- 
served transition discs in regards to gap sizes and mass ac- 
cretion rates. However, our model shows disc masses (we 
take masses outside gaps) an order of magnitude smaller 
than those for the most massive observed transition discs. 
This may indicate that the dust-to-gas ratios are large in 
the outer discs of observed transition discs while our model 
assumes a fixed dust opacity. 
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